Bernoulli Map
   HOME

TheInfoList



OR:

The dyadic transformation (also known as the dyadic map, bit shift map, 2''x'' mod 1 map, Bernoulli map, doubling map or sawtooth map) is the mapping (i.e.,
recurrence relation In mathematics, a recurrence relation is an equation according to which the nth term of a sequence of numbers is equal to some combination of the previous terms. Often, only k previous terms of the sequence appear in the equation, for a parameter ...
) : T: , 1) \to [0, 1)^\infty : x \mapsto (x_0, x_1, x_2, \ldots) (where [0, 1)^\infty is the set of sequences from [0, 1)) produced by the rule : x_0 = x : \text n \ge 0,\ x_ = (2 x_n) \bmod 1. Equivalently, the dyadic transformation can also be defined as the iterated function map of the piecewise linear function : T(x)=\begin2x & 0 \le x < \frac \\2x-1 & \frac \le x < 1. \end The name ''bit shift map'' arises because, if the value of an iterate is written in
binary Binary may refer to: Science and technology Mathematics * Binary number, a representation of numbers using only two digits (0 and 1) * Binary function, a function that takes two arguments * Binary operation, a mathematical operation that ta ...
notation, the next iterate is obtained by shifting the binary point one bit to the right, and if the bit to the left of the new binary point is a "one", replacing it with a zero. The dyadic transformation provides an example of how a simple 1-dimensional map can give rise to
chaos Chaos or CHAOS may refer to: Arts, entertainment and media Fictional elements * Chaos (''Kinnikuman'') * Chaos (''Sailor Moon'') * Chaos (''Sesame Park'') * Chaos (''Warhammer'') * Chaos, in ''Fabula Nova Crystallis Final Fantasy'' * Cha ...
. This map readily generalizes to several others. An important one is the
beta transformation Beta (, ; uppercase , lowercase , or cursive Greek, cursive ; grc, βῆτα, bē̂ta or ell, βήτα, víta) is the second letter of the Greek alphabet. In the system of Greek numerals, it has a value of 2. In Modern Greek, it represents th ...
, defined as T_\beta (x)=\beta x\bmod 1. This map has been extensively studied by many authors. It was introduced by Alfréd Rényi in 1957, and an invariant measure for it was given by
Alexander Gelfond Alexander Osipovich Gelfond (russian: Алекса́ндр О́сипович Ге́льфонд; 24 October 1906 – 7 November 1968) was a Soviet mathematician. Gelfond's theorem, also known as the Gelfond-Schneider theorem is named after hi ...
in 1959 and again independently by Bill Parry in 1960.


Relation to the Bernoulli process

The map can be obtained as a
homomorphism In algebra, a homomorphism is a structure-preserving map between two algebraic structures of the same type (such as two groups, two rings, or two vector spaces). The word ''homomorphism'' comes from the Ancient Greek language: () meaning "same" ...
on the Bernoulli process. Let \Omega = \^ be the set of all semi-infinite strings of the letters H and T. These can be understood to be the flips of a coin, coming up heads or tails. Equivalently, one can write \Omega = \^ the space of all (semi-)infinite strings of binary bits. The word "infinite" is qualified with "semi-", as one can also define a different space \^ consisting of all doubly-infinite (double-ended) strings; this will lead to the Baker's map. The qualification "semi-" is dropped below. This space has a natural shift operation, given by :T(b_0, b_1, b_2, \dots) = (b_1, b_2, \dots) where (b_0, b_1, \dots) is an infinite string of binary digits. Given such a string, write :x = \sum_^\infty \frac. The resulting x is a
real number In mathematics, a real number is a number that can be used to measure a ''continuous'' one-dimensional quantity such as a distance, duration or temperature. Here, ''continuous'' means that values can have arbitrarily small variations. Every ...
in the
unit interval In mathematics, the unit interval is the closed interval , that is, the set of all real numbers that are greater than or equal to 0 and less than or equal to 1. It is often denoted ' (capital letter ). In addition to its role in real analysis ...
0 \le x \le 1. The shift T induces a
homomorphism In algebra, a homomorphism is a structure-preserving map between two algebraic structures of the same type (such as two groups, two rings, or two vector spaces). The word ''homomorphism'' comes from the Ancient Greek language: () meaning "same" ...
, also called T, on the unit interval. Since T(b_0, b_1, b_2, \dots) = (b_1, b_2, \dots), one can easily see that T(x)=2x\bmod 1. For the doubly-infinite sequence of bits \Omega = 2^, the induced homomorphism is the Baker's map. The dyadic sequence is then just the sequence :(x, T(x), T^2(x), T^3(x), \dots) That is, x_n = T^n(x).


The Cantor set

Note that the sum :y=\sum_^\infty \frac gives the
Cantor function In mathematics, the Cantor function is an example of a function that is continuous, but not absolutely continuous. It is a notorious counterexample in analysis, because it challenges naive intuitions about continuity, derivative, and measure. ...
, as conventionally defined. This is one reason why the set \^\mathbb is sometimes called the Cantor set.


Rate of information loss and sensitive dependence on initial conditions

One hallmark of chaotic dynamics is the loss of information as simulation occurs. If we start with information on the first ''s'' bits of the initial iterate, then after ''m'' simulated iterations (''m'' < ''s'') we only have ''s'' − ''m'' bits of information remaining. Thus we lose information at the exponential rate of one bit per iteration. After ''s'' iterations, our simulation has reached the fixed point zero, regardless of the true iterate values; thus we have suffered a complete loss of information. This illustrates sensitive dependence on initial conditions—the mapping from the truncated initial condition has deviated exponentially from the mapping from the true initial condition. And since our simulation has reached a fixed point, for almost all initial conditions it will not describe the dynamics in the qualitatively correct way as chaotic. Equivalent to the concept of information loss is the concept of information gain. In practice some real-world process may generate a sequence of values (''x''''n'') over time, but we may only be able to observe these values in truncated form. Suppose for example that ''x''0 = 0.1001101, but we only observe the truncated value 0.1001. Our prediction for ''x''1 is 0.001. If we wait until the real-world process has generated the true ''x''1 value 0.001101, we will be able to observe the truncated value 0.0011, which is more accurate than our predicted value 0.001. So we have received an information gain of one bit.


Relation to tent map and logistic map

The dyadic transformation is
topologically semi-conjugate In mathematics, two functions are said to be topologically conjugate if there exists a homeomorphism that will conjugate the one into the other. Topological conjugacy, and related-but-distinct of flows, are important in the study of iterated fun ...
to the unit-height
tent map A tent () is a shelter consisting of sheets of fabric or other material draped over, attached to a frame of poles or a supporting rope. While smaller tents may be free-standing or attached to the ground, large tents are usually anchored using g ...
. Recall that the unit-height tent map is given by :x_ = f_1(x_n) = \begin x_n & \mathrm~~ x_n \le 1/2 \\ 1-x_n & \mathrm~~ x_n \ge 1/2 \end The conjugacy is explicitly given by :S(x)=\sin \pi x so that :f_1 = S^ \circ T \circ S That is, f_1(x) = S^(T(S(x))). This is stable under iteration, as :f_1^n = f_1\circ\cdots\circ f_1 = S^ \circ T \circ S \circ S^ \circ \cdots \circ T \circ S = S^ \circ T^n \circ S It is also conjugate to the chaotic ''r'' = 4 case of the logistic map. The ''r'' = 4 case of the logistic map is z_=4z_n(1-z_n); this is related to the
bit shift In computer programming, a bitwise operation operates on a bit string, a bit array or a binary numeral (considered as a bit string) at the level of its individual bits. It is a fast and simple action, basic to the higher-level arithmetic oper ...
map in variable ''x'' by :z_n =\sin^2 (2 \pi x_n). There is also a semi-conjugacy between the dyadic transformation (here named angle doubling map) and the
quadratic polynomial In mathematics, a quadratic polynomial is a polynomial of degree two in one or more variables. A quadratic function is the polynomial function defined by a quadratic polynomial. Before 20th century, the distinction was unclear between a polynomial ...
. Here, the map doubles angles measured in turns. That is, the map is given by :\theta\mapsto 2\theta\bmod 2\pi.


Periodicity and non-periodicity

Because of the simple nature of the dynamics when the iterates are viewed in binary notation, it is easy to categorize the dynamics based on the initial condition: If the initial condition is
irrational Irrationality is cognition, thinking, talking, or acting without inclusion of rationality. It is more specifically described as an action or opinion given through inadequate use of reason, or through emotional distress or cognitive deficiency. T ...
(as almost all points in the unit interval are), then the dynamics are non-periodic—this follows directly from the definition of an irrational number as one with a non-repeating binary expansion. This is the chaotic case. If ''x''0 is
rational Rationality is the quality of being guided by or based on reasons. In this regard, a person acts rationally if they have a good reason for what they do or a belief is rational if it is based on strong evidence. This quality can apply to an abi ...
the image of ''x''0 contains a finite number of distinct values within forward_orbit_of_''x''0_is_eventually_periodic,_with_period_equal_to_the_period_of_the_
forward_orbit_of_''x''0_is_eventually_periodic,_with_period_equal_to_the_period_of_the_Binary_numeral_system">binary_ Binary_may_refer_to: __Science_and_technology___Mathematics_ *_Binary_number,_a_representation_of_numbers_using_only_two_digits_(0_and_1) *_Binary_function,_a_function_that_takes_two_arguments *_Binary_operation,_a_mathematical_operation_that_ta_...
_expansion_of_''x''0.__Specifically,_if_the_initial_condition_is_a_rational_number_with_a_finite_binary_expansion_of_''k''_bits,_then_after_''k''_iterations_the_iterates_reach_the_fixed_point_0; if_the_initial_condition_is_a_rational_number_with_a_''k''-bit_transient_(''k'' ≥ 0)_followed_by_a_''q''-bit_sequence_(''q'' > 1)_that_repeats_itself_infinitely,_then_after_''k''_iterations_the_iterates_reach_a_cycle_of_length ''q''._Thus_cycles_of_all_lengths_are_possible. For_example,_the_forward_orbit_of_11/24_is: :_\frac_\mapsto_\frac_\mapsto_\frac_\mapsto_\frac_\mapsto_\frac_\mapsto_\frac_\mapsto_\frac_\mapsto_\cdots,_ which_has_reached_a_cycle_of_period_2.__Within_any_subinterval_of_[0, 1),_no_matter_how_small,_there_are_therefore_an_infinite_number_of_points_whose_orbits_are_eventually_periodic,_and_an_infinite_number_of_points_whose_orbits_are_never_periodic._This_sensitive_dependence_on_initial_conditions_is_a_characteristic_of_list_of_chaotic_maps.html" "title="Binary_numeral_system.html" "title="orbit_(dynamics).html" ;"title=", 1) and the orbit (dynamics)">forward orbit of ''x''0 is eventually periodic, with period equal to the period of the Binary numeral system">binary Binary may refer to: Science and technology Mathematics * Binary number, a representation of numbers using only two digits (0 and 1) * Binary function, a function that takes two arguments * Binary operation, a mathematical operation that ta ...
expansion of ''x''0. Specifically, if the initial condition is a rational number with a finite binary expansion of ''k'' bits, then after ''k'' iterations the iterates reach the fixed point 0; if the initial condition is a rational number with a ''k''-bit transient (''k'' ≥ 0) followed by a ''q''-bit sequence (''q'' > 1) that repeats itself infinitely, then after ''k'' iterations the iterates reach a cycle of length ''q''. Thus cycles of all lengths are possible. For example, the forward orbit of 11/24 is: : \frac \mapsto \frac \mapsto \frac \mapsto \frac \mapsto \frac \mapsto \frac \mapsto \frac \mapsto \cdots, which has reached a cycle of period 2. Within any subinterval of [0, 1), no matter how small, there are therefore an infinite number of points whose orbits are eventually periodic, and an infinite number of points whose orbits are never periodic. This sensitive dependence on initial conditions is a characteristic of list of chaotic maps">chaotic maps Chaotic was originally a Danish trading card game. It expanded to an online game in America which then became a television program based on the game. The program was able to be seen on 4Kids TV (Fox affiliates, nationwide), Jetix, The CW4Kid ...
.


Periodicity via bit shifts

The periodic and non-periodic orbits can be more easily understood not by working with the map T(x)=2x\bmod 1 directly, but rather with the
bit shift In computer programming, a bitwise operation operates on a bit string, a bit array or a binary numeral (considered as a bit string) at the level of its individual bits. It is a fast and simple action, basic to the higher-level arithmetic oper ...
map T(b_0,b_1,b_2,\dots) = (b_1, b_2,\dots) defined on the
Cantor space In mathematics, a Cantor space, named for Georg Cantor, is a topological abstraction of the classical Cantor set: a topological space is a Cantor space if it is homeomorphic to the Cantor set. In set theory, the topological space 2ω is called "the ...
\Omega=\^\mathbb. That is, the
homomorphism In algebra, a homomorphism is a structure-preserving map between two algebraic structures of the same type (such as two groups, two rings, or two vector spaces). The word ''homomorphism'' comes from the Ancient Greek language: () meaning "same" ...
:x=\sum_^\infty \frac is basically a statement that the Cantor set can be mapped into the reals. It is a
surjection In mathematics, a surjective function (also known as surjection, or onto function) is a function that every element can be mapped from element so that . In other words, every element of the function's codomain is the image of one element of ...
: every dyadic rational has not one, but two distinct representations in the Cantor set. For example, :0.1000000\dots = 0.011111\dots This is just the binary-string version of the famous 0.999... = 1 problem. The doubled representations hold in general: for any given finite-length initial sequence b_0,b_1,b_2,\dots,b_ of length k, one has :b_0,b_1,b_2,\dots,b_,1,0,0,0,\dots = b_0,b_1,b_2,\dots,b_,0,1,1,1,\dots The initial sequence b_0,b_1,b_2,\dots,b_ corresponds to the non-periodic part of the orbit, after which iteration settles down to all zeros (equivalently, all-ones). Expressed as bit strings, the periodic orbits of the map can be seen to the rationals. That is, after an initial "chaotic" sequence of b_0,b_1,b_2,\dots,b_, a periodic orbit settles down into a repeating string b_k,b_,b_,\dots,b_ of length m. It is not hard to see that such repeating sequences correspond to rational numbers. Writing :y = \sum_^ b_2^ one then clearly has :\sum_^\infty b_2^ = y\sum_^\infty 2^ = \frac Tacking on the initial non-repeating sequence, one clearly has a rational number. In fact, ''every'' rational number can be expressed in this way: an initial "random" sequence, followed by a cycling repeat. That is, the periodic orbits of the map are in one-to-one correspondence with the rationals. This phenomenon is note-worthy, because something similar happens in many chaotic systems. For example, geodesics on
compact Compact as used in politics may refer broadly to a pact or treaty; in more specific cases it may refer to: * Interstate compact * Blood compact, an ancient ritual of the Philippines * Compact government, a type of colonial rule utilized in British ...
manifolds can have periodic orbits that behave in this way. Keep in mind, however, that the rationals are a set of
measure zero In mathematical analysis, a null set N \subset \mathbb is a measurable set that has measure zero. This can be characterized as a set that can be covered by a countable union of intervals of arbitrarily small total length. The notion of null ...
in the reals. Almost all orbits are ''not'' periodic! The aperiodic orbits correspond to the irrational numbers. This property also holds true in a more general setting. An open question is to what degree the behavior of the periodic orbits constrain the behavior of the system as a whole. Phenomena such as
Arnold diffusion In applied mathematics, Arnold diffusion is the phenomenon of instability of integrable Hamiltonian systems. The phenomenon is named after Vladimir Arnold who was the first to publish a result in the field in 1964. More precisely, Arnold diffusio ...
suggest that the general answer is "not very much".


Density formulation

Instead of looking at the orbits of individual points under the action of the map, it is equally worthwhile to explore how the map affects densities on the unit interval. That is, imagine sprinkling some dust on the unit interval; it is denser in some places than in others. What happens to this density as one iterates? Write \rho: ,1to\mathbb as this density, so that x\mapsto\rho(x). To obtain the action of T on this density, one needs to find all points y=T^(x) and write Dean J. Driebe, Fully Chaotic Maps and Broken Time Symmetry, (1999) Kluwer Academic Publishers, Dordrecht Netherlands :\rho(x) \mapsto \sum_ \frac The denominator in the above is the
Jacobian determinant In vector calculus, the Jacobian matrix (, ) of a vector-valued function of several variables is the matrix of all its first-order partial derivatives. When this matrix is square, that is, when the function takes the same number of variables ...
of the transformation, here it is just the
derivative In mathematics, the derivative of a function of a real variable measures the sensitivity to change of the function value (output value) with respect to a change in its argument (input value). Derivatives are a fundamental tool of calculus. ...
of T and so T^\prime(y)=2. Also, there are obviously only two points in the preimage of T^(x), these are y=x/2 and y=(x+1)/2. Putting it all together, one gets :\rho(x) \mapsto \frac\rho\!\left(\frac\right) + \frac\rho\!\left(\frac\right) By convention, such maps are denoted by \mathcal so that in this case, write :\left mathcal _T\rho\rightx) = \frac\rho\!\left(\frac\right) + \frac\rho\!\left(\frac\right) The map \mathcal_T is a linear operator, as one easily sees that \mathcal_T(f+g)= \mathcal_T(f) + \mathcal_T(g) and \mathcal_T(af)= a\mathcal_T(f) for all functions f,g on the unit interval, and all constants a. Viewed as a linear operator, the most obvious and pressing question is: what is its
spectrum A spectrum (plural ''spectra'' or ''spectrums'') is a condition that is not limited to a specific set of values but can vary, without gaps, across a continuum. The word was first used scientifically in optics to describe the rainbow of colors ...
? One
eigenvalue In linear algebra, an eigenvector () or characteristic vector of a linear transformation is a nonzero vector that changes at most by a scalar factor when that linear transformation is applied to it. The corresponding eigenvalue, often denoted ...
is obvious: if \rho(x)=1 for all x then one obviously has \mathcal_T\rho=\rho so the uniform density is invariant under the transformation. This is in fact the largest eigenvalue of the operator \mathcal_T, it is the Frobenius–Perron eigenvalue. The uniform density is, in fact, nothing other than the
invariant measure In mathematics, an invariant measure is a measure that is preserved by some function. The function may be a geometric transformation. For examples, circular angle is invariant under rotation, hyperbolic angle is invariant under squeeze mapping, ...
of the dyadic transformation. To explore the spectrum of \mathcal_T in greater detail, one must first limit oneself to a suitable
space of functions In mathematics, a function space is a set of functions between two fixed sets. Often, the domain and/or codomain will have additional structure which is inherited by the function space. For example, the set of functions from any set into a ve ...
(on the unit interval) to work with. This might be the space of Lebesgue measurable functions, or perhaps the space of
square integrable In mathematics, a square-integrable function, also called a quadratically integrable function or L^2 function or square-summable function, is a real- or complex-valued measurable function for which the integral of the square of the absolute value i ...
functions, or perhaps even just
polynomial In mathematics, a polynomial is an expression consisting of indeterminates (also called variables) and coefficients, that involves only the operations of addition, subtraction, multiplication, and positive-integer powers of variables. An example ...
s. Working with any of these spaces is surprisingly difficult, although a spectrum can be obtained.


Borel space

A vast amount of simplification results if one instead works with the
Cantor space In mathematics, a Cantor space, named for Georg Cantor, is a topological abstraction of the classical Cantor set: a topological space is a Cantor space if it is homeomorphic to the Cantor set. In set theory, the topological space 2ω is called "the ...
\Omega=\^\mathbb, and functions \rho:\Omega\to\mathbb. Some caution is advised, as the map T(x)=2x\bmod 1 is defined on the
unit interval In mathematics, the unit interval is the closed interval , that is, the set of all real numbers that are greater than or equal to 0 and less than or equal to 1. It is often denoted ' (capital letter ). In addition to its role in real analysis ...
of the
real number line In elementary mathematics, a number line is a picture of a graduated straight line that serves as visual representation of the real numbers. Every point of a number line is assumed to correspond to a real number, and every real number to a poin ...
, assuming the
natural topology In any domain of mathematics, a space has a natural topology if there is a topology on the space which is "best adapted" to its study within the domain in question. In many cases this imprecise definition means little more than the assertion that ...
on the reals. By contrast, the map T(b_0, b_1, b_2, \dots)=(b_1, b_2, \dots) is defined on the
Cantor space In mathematics, a Cantor space, named for Georg Cantor, is a topological abstraction of the classical Cantor set: a topological space is a Cantor space if it is homeomorphic to the Cantor set. In set theory, the topological space 2ω is called "the ...
\Omega = \^, which by convention is given a very different
topology In mathematics, topology (from the Greek words , and ) is concerned with the properties of a geometric object that are preserved under continuous deformations, such as stretching, twisting, crumpling, and bending; that is, without closing ...
, the
product topology In topology and related areas of mathematics, a product space is the Cartesian product of a family of topological spaces equipped with a natural topology called the product topology. This topology differs from another, perhaps more natural-s ...
. There is a potential clash of topologies; some care must be taken. However, as presented above, there is a homomorphism from the Cantor set into the reals; fortunately, it maps
open set In mathematics, open sets are a generalization of open intervals in the real line. In a metric space (a set along with a distance defined between any two points), open sets are the sets that, with every point , contain all points that are su ...
s into open sets, and thus preserves notions of continuity. To work with the Cantor set \Omega=\^, one must provide a topology for it; by convention, this is the
product topology In topology and related areas of mathematics, a product space is the Cartesian product of a family of topological spaces equipped with a natural topology called the product topology. This topology differs from another, perhaps more natural-s ...
. By adjoining set-complements, it can be extended to a Borel space, that is, a sigma algebra. The topology is that of cylinder sets. A cylinder set has the generic form :(*,*,*,\dots,*,b_k,b_,*,\dots, *,b_m,*,\dots) where the * are arbitrary bit values (not necessarily all the same), and the b_k, b_m, \dots are a finite number of specific bit-values scattered in the infinite bit-string. These are the open sets of the topology. The canonical measure on this space is the Bernoulli measure for the fair coin-toss. If there is just one bit specified in the string of arbitrary positions, the measure is 1/2. If there are two bits specified, the measure is 1/4, and so on. One can get fancier: given a real number 0 < p < 1 one can define a measure :\mu_p( *,\dots,*,b_k,*,\dots) = p^n(1-p)^m if there are n heads and m tails in the sequence. The measure with p=1/2 is preferred, since it is preserved by the map :(b_0, b_1, b_2, \dots) \mapsto x = \sum_^\infty \frac. So, for example, (0,*,\cdots) maps to the interval ,1/2/math> and (1,*,\dots) maps to the interval /2,1/math> and both of these intervals have a measure of 1/2. Similarly, (*,0,*,\dots) maps to the interval ,1/4cup /2,3/4/math> which still has the measure 1/2. That is, the embedding above preserves the measure. An alternative is to write :(b_0, b_1, b_2, \dots) \mapsto x = \sum_^\infty \left _n p^ + (1-b_n)(1-p)^\right/math> which preserves the measure \mu_p. That is, it maps such that the measure on the unit interval is again the Lebesgue measure.


Frobenius–Perron operator

Denote the collection of all open sets on the Cantor set by \mathcal and consider the set \mathcal of all arbitrary functions f:\mathcal\to\mathbb. The shift T induces a
pushforward The notion of pushforward in mathematics is "dual" to the notion of pullback, and can mean a number of different but closely related things. * Pushforward (differential), the differential of a smooth map between manifolds, and the "pushforward" op ...
:f\circ T^ defined by \left(f \circ T^\right)\!(x) = f(T^(x)). This is again some function \mathcal\to\mathbb. In this way, the map T induces another map \mathcal_T on the space of all functions \mathcal\to\mathbb. That is, given some f:\mathcal\to\mathbb, one defines :\mathcal_T f = f \circ T^ This linear operator is called the
transfer operator Transfer may refer to: Arts and media * ''Transfer'' (2010 film), a German science-fiction movie directed by Damir Lukacevic and starring Zana Marjanović * ''Transfer'' (1966 film), a short film * ''Transfer'' (journal), in management studies ...
or the ''Ruelle–Frobenius–Perron operator''. The largest eigenvalue is the Frobenius–Perron eigenvalue, and in this case, it is 1. The associated
eigenvector In linear algebra, an eigenvector () or characteristic vector of a linear transformation is a nonzero vector that changes at most by a scalar factor when that linear transformation is applied to it. The corresponding eigenvalue, often denoted ...
is the invariant measure: in this case, it is the Bernoulli measure. Again, \mathcal_T(\rho)= \rho when \rho(x)=1.


Spectrum

To obtain the spectrum of \mathcal_T, one must provide a suitable set of basis functions for the space \mathcal. One such choice is to restrict \mathcal to the set of all polynomials. In this case, the operator has a
discrete spectrum A observable, physical quantity is said to have a discrete spectrum if it takes only distinct values, with gaps between one value and the next. The classical example of discrete spectrum (for which the term was first used) is the characterist ...
, and the
eigenfunction In mathematics, an eigenfunction of a linear operator ''D'' defined on some function space is any non-zero function f in that space that, when acted upon by ''D'', is only multiplied by some scaling factor called an eigenvalue. As an equation, th ...
s are (curiously) the
Bernoulli polynomial In mathematics, the Bernoulli polynomials, named after Jacob Bernoulli, combine the Bernoulli numbers and binomial coefficients. They are used for series expansion of functions, and with the Euler–MacLaurin formula. These polynomials occur ...
s! (This coincidence of naming was presumably not known to Bernoulli.) Indeed, one can easily verify that :\mathcal_T B_n= 2^B_n where the B_n are the Bernoulli polynomials. This follows because the Bernoulli polynomials obey the identity :\fracB_n\!\left(\frac\right) + \fracB_n\!\left(\frac\right) = 2^B_n(y) Note that B_0(x)=1. Another basis is provided by the Haar basis, and the functions spanning the space are the
Haar wavelet In mathematics, the Haar wavelet is a sequence of rescaled "square-shaped" functions which together form a wavelet family or basis. Wavelet analysis is similar to Fourier analysis in that it allows a target function over an interval to be repres ...
s. In this case, one finds a continuous spectrum, consisting of the unit disk on the complex plane. Given z\in\mathbb in the unit disk, so that , z, <1, the functions :\psi_(x)=\sum_^\infty z^n \exp i\pi(2k+1)2^nx obey :\mathcal_T \psi_= z\psi_ for k\in\mathbb. This is a complete basis, in that every
integer An integer is the number zero (), a positive natural number (, , , etc.) or a negative integer with a minus sign ( −1, −2, −3, etc.). The negative numbers are the additive inverses of the corresponding positive numbers. In the languag ...
can be written in the form (2k+1)2^n. The Bernoulli polynomials are recovered by setting k=0 and z=\frac, \frac, \dots A complete basis can be given in other ways, as well; they may be written in terms of the
Hurwitz zeta function In mathematics, the Hurwitz zeta function is one of the many zeta functions. It is formally defined for complex variables with and by :\zeta(s,a) = \sum_^\infty \frac. This series is absolutely convergent for the given values of and and c ...
. Another complete basis is provided by the Takagi function. This is a fractal, differentiable-nowhere function. The eigenfunctions are explicitly of the form :\mbox_(x) = \sum_^\infty w^n s((2k+1)2^x) where s(x) is the
triangle wave A triangular wave or triangle wave is a non-sinusoidal waveform named for its triangular shape. It is a periodic, piecewise linear, continuous real function. Like a square wave, the triangle wave contains only odd harmonics. However, ...
. One has, again, :\mathcal_T \mbox_ = w\;\mbox_. All of these different bases can be expressed as linear combinations of one-another. In this sense, they are equivalent. The fractal eigenfunctions show an explicit symmetry under the fractal
groupoid In mathematics, especially in category theory and homotopy theory, a groupoid (less often Brandt groupoid or virtual group) generalises the notion of group in several equivalent ways. A groupoid can be seen as a: *'' Group'' with a partial func ...
of the modular group; this is developed in greater detail in the article on the Takagi function (the blancmange curve). Perhaps not a surprise; the Cantor set has exactly the same set of symmetries (as do the
continued fraction In mathematics, a continued fraction is an expression obtained through an iterative process of representing a number as the sum of its integer part and the reciprocal of another number, then writing this other number as the sum of its integer ...
s.) This then leads elegantly into the theory of
elliptic equation An elliptic equation can mean: * The equation of an ellipse * An elliptic curve, describing the relationships between invariants of an ellipse * A differential equation with an elliptic operator * An elliptic partial differential equation Second ...
s and modular forms.


Relation to the Ising model

The Hamiltonian of the zero-field one-dimensional
Ising model The Ising model () (or Lenz-Ising model or Ising-Lenz model), named after the physicists Ernst Ising and Wilhelm Lenz, is a mathematical model of ferromagnetism in statistical mechanics. The model consists of discrete variables that represent ...
of 2N spins with periodic boundary conditions can be written as :H(\sigma) = g \sum_\sigma_i\sigma_. Letting C be a suitably chosen normalization constant and \beta be the inverse temperature for the system, the partition function for this model is given by :Z = \sum_\prod_Ce^. We can implement the
renormalization group In theoretical physics, the term renormalization group (RG) refers to a formal apparatus that allows systematic investigation of the changes of a physical system as viewed at different scales. In particle physics, it reflects the changes in the ...
by integrating out every other spin. In so doing, one finds that Z can also be equated with the partition function for a smaller system with but N spins, :Z = \sum_\prod_\mathcal ^, provided we replace C and \beta g with renormalized values \mathcal /math> and \mathcal beta g/math> satisfying the equations :\mathcal 2= 4\cosh(2\beta g)C^4, :e^= \cosh(2\beta g). Suppose now that we allow \beta g to be complex and that \operatorname \beta g\frac+\pi n for some n\in \mathbb. In that case we can introduce a parameter t\in[0, 1) related to \beta g via the equation :e^= i\tan\big(\pi(t-\frac)\big), and the resulting renormalization group transformation for t will be precisely the dyadic map: M. Bosschaert; C. Jepsen; F. Popov, “Chaotic RG flow in tensor models”, Physical Review D, 105, 2022, p. 065021. :\mathcal[t]=2t \bmod 1 .


See also

* Bernoulli process * Bernoulli scheme * Gilbert–Shannon–Reeds model, a random distribution on permutations given by applying the doubling map to a set of ''n'' uniformly random points on the unit interval


Notes


References

* Dean J. Driebe, ''Fully Chaotic Maps and Broken Time Symmetry'', (1999) Kluwer Academic Publishers, Dordrecht Netherlands * Linas Vepstas,
The Bernoulli Map, the Gauss-Kuzmin-Wirsing Operator and the Riemann Zeta
', (2004) {{Chaos theory Chaotic maps